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Abstract 

In spite of its long history and classical character which goes back even 
to d'Alembert and Lagrange, the problems of constraints in mechanics of 
continua is still mysterious and full of misunderstandings. Let us mention 
the problem of difference between special solutions of unconstrained systems 
and constraints dynamics, holonomic and nonholonomic constrained, various 
versions of so-called vakonomy, sub-Riemannian structures, servoconstraints, 
and programme motion. There are some strange confusions between field the- 
ory and mechanics of continua; many problems call for conceptual cleaning, 
although both the mechanics of continua and field theory are formally based 
on partial differential equations. The author is deeply indebted to professor 
Czeslaw Wozniak for his help during many years, and deep understanding 
and support 

1 Constrained dynamics versus constrained solu- 
tions 

The problem of constraints in mechanics is very delicate and full of confusions 
and misconceptions. Situation is maximally clear in mechanics of system with 
a finite number of degrees of freedom, subject to holonomic constraints. Let us 
consider a system of N material points moving in n-dimensional Euclidean space 
M. "Physically" n — 3, but let us forget about this. If some other symbols turn 
out to be needed, the ./V-dimensional linear space of translations in M will be 
denoted by V , and the metric tensor by g £ V* ® V*; obviously, g is symmetric 
and positively definite. 

Let the force acting on the A-th material point be denoted by 

F A (Ti, . . . ,7n;vi, . . . ,v N ;t) 
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the explicit dependence on time t is admissible. Obviously, r a is the radius-vector 
of the A-th material point (with respect to some fixed origin in M) and va — dfA/dt 
is its velocity vector (independent on the origin, but evidently dependent on the 
choice of the inertial reference frame) . Unconstrained equations of motion have the 
obvious Newton form: 



d r a -pr / _ _ dTi dr 



Kinetic energy is given by 



1 ^ dr A df A , n , 
a ra A T ' T (2) 



T ..... 

2 ^ dt dt 

A 

with the obvious symbol for the scalar product: 

u-v = g(u,v) — gijU l w J . (3) 

If we use generalized coordinates Q a and the corresponding kinematical metric 
tensor G on the configuration space Q :— M N , 

dr A dfA ,.s 

and its contravariant inverse G ab , 

G ac G cb = 6\, (5) 

equations of motion become 

^(..., v,.44 

Here generalized accelerations are given by 

dp ~ dt 2 + bc dt dt ' [ ' 

r a f, c are Christoffcl coefficients for G, 

T a fcc = 2^* al (Gib.c + Gi c ^b — Gbc,i) , (8) 

comma denotes the partial differentiation with respect to Q a , and generalized forces 
are given by 

A 

Obviously, the configuration kinematical metric G is here flat, i.e., its curvature 
tensor does vanish. Nevertheless, equations ([6]), (0 are more general and hold also 
for the general Riemann space (Q, G), used as a configuration space. 
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For the potential system, the forces Fa are given by 



dV ■ dV 

F 1 a = -9 n tt— = -g v jfj-5 B A. (10) 
or J a or -3 b 



In terms of generalized coordinates: 

u dV 

pa =- G W 

Equations of motion are then derivable from the variational principle: 

61 = 6 [ Ldt = 0, (12) 



0; (13) 



where L = T — V, and have the Euler-Lagrange form 

SI dL D dL 



the symbol D/Dt denoting the total (substantial) derivative with respect to time 
t; do not confuse it with the covariant differentiation. 

If in addition to the Lagrangian background some other forces 4> a , e.g., dissipa- 
tive ones (friction) are present, then (TlB")) is replaced by 

(jab? ■ (14) 



dQ a Dt dQ a 

The problem of constraints is as follows: In addition to explicitly given forces 
like Fa, F a , etc., there are additional ones, which confine motion to some sub- 
manifold W C Q of dimension / < dimQ. More precisely, this is the problem 
of holonomic constraints. The manifold W may be analytically described by the 
system of equations: 

F a (...,Q b ,...) = 0, o = l,...,m = dimQ-/. (15) 

And now some very important point comes, namely one concerning the distinc- 
tion between two different problems: 

Procedure 1. Problem of special solutions of (|13l) / (fT4|) satisfying equations 
(fT5|) , i.e., placed on 

W :={zeQ:F a {z) = 0,a = l,...,m} (16) 

(it is implicitly assumed that <p a are functionally independent, at least in a neigh- 
bourhood of W). 

Procedure 2. Problem (|T3"f / (fl4|) dynamically modified by (|15p . This problem 
consists in that in addition to a priori given forces (jTJ), (J9j), (flO|) . etc., there are some 
additional ones, usually even not described explicitly (usually, but not always), 
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TZa, TZ a which result in small-amplitude (but sometimes large-energy and high- 
frequency) oscillations vertical to W. Their orthogonal projections onto W do not 
satisfy original equations, but something else. 

Obviously, it is well-known that the total, joint system consisting of the pri- 
mary equations of motion in any of the versions (JXJ) , (|13[) . (|14p and constraining 
conditions (|15[) is inconsistent: as a rule, it is intrinsically contradictory. There 
are only exceptional situations when inconsistency does not appear. For example, 
when the forces Fa do vanish, i.e., motion is free, and W is an affine subman- 
ifold of the Euclidean space Q (functions F a are affine, i.e., roughly speaking, 
linear- inhomogeneous), then all straight-lines swept with constant velocities are 
solutions (they form the general solution of the problem of motion in W). But 
of course, the situation changes drastically when the general constraints of rigid 
motion are imposed onto the free motion problem. Then the only solutions of the 
above Procedure 1 are the total free translations. But we are aware that with a 
good approximation the general rigid motion with nontrivial rotations does exist 
and is both theoretically interesting and practically relevant; one knows that from 
elementary school. The point is that one deals then with the above Procedure 
2. The solution is classical, known from ages and belongs mainly to Lagrange 
and d'Alembert. Nevertheless, there are situations when this classical procedure is 
forgotten. From some point of view this is the problem of passive control, based 
on the natural contact-rolling-friction mechanism, without servomechanism-aided 
and computer-aided procedures and active control elements, cf PQ-[B]. From some 
point of view the idea is peculiar. Namely, instead of real auxiliary forces or- 
thogonal to W, defined around W, and just vanishing at W itself, one introduced 
symbolic d'Alembert reactions which are just defined on W itself and maintain the 
among-constraints motion. Therefore, ([!]) is replaced by 

d 2 TA ( _ _ dri dr N \ — (_ _ dri rffjv . 

mA ^ = FA { ri '-- rN; ^r--^r ;t ) +1ZA ( v ri '---' rjv; ^'---^r ;i 

(17) 

and this equation is completed and treated jointly with constraint equations (|15|) 

F a (r 1 ,...,r N )=F a (...,Q b ,...) = (18) 

But of course the joint system (TT71) k (IT8l) is over-determined unless some con- 
stitutive conditions are imposed on TZa- The d'Alembert-Lagrange procedure of 
passive control tells us that the along-constraints motion is neither affected nor 
maintained without the energy transfer to the system, i.e., reactions TZa do not 
do any work, i.e., their mechanical power does vanish; they are orthogonal to the 
constraints surface W, i.e., to all virtual motions: 

^ = £7^-^ = (19) 

A 

for any functions R 9 t —> r a (t) satisfying equations: 

F a (ri (t),...,r N (t)) = identically over t, (20) 
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therefore, for any virtual velocities subject to 

Obviously, the symbols dF a /dfA are abbreviations for the systems of derivatives 
with respect to the components r % a- This means that reactions TZa are given by: 

^ = Y>'w a (22) 

a 

with the same meaning of dF a /dfA as above. Here A a are some apriori unknown 
Lagrange multipliers. The total system of equations of motion in the sense of the 
Procedure 2 above, consists of (|17l) . ([20]) . (|2"2"j) . This is the system of / + m = 
dim Q+m equations for the system of dim Q+m variables 7a, A a . Eliminating A a — s 
we obtain the system of / independent equations for / independent generalized 
coordinates parametrizing W C Q. The knowledge of A a , i.e., reactions TZa is 
essential for the design of the system of factors maintaining the constraints. 

Let us go back to a more homogeneous notation. We consider analytical me- 
chanics in a differential manifold Q of dimension dimQ = k, with the Lagrangian 
background L : TQ — >• K ( TQ denoting the tangent bundle of Q, i.e., manifold of 
generalized coordinates q l and velocities v l = <?'), and with certain non-Lagrangian 
forces, e.g., dissipative ones, represented in the covariant terms as Di. If Q is en- 
dowed with Riemannian structure G, and something like the "magnetic field" Ai 
(covector potential) is present, and the scalar potential V as well, then the typical 
school example is: 

L (q, v) = \G l3 (q) wV + eAi (q) v* + V (q) . (23) 

Typical model of dissipative forces is 

D i (q,v) = -d ij (q,v)v j , (24) 

where dij (q, 0) = 0, and dij is symmetric and positively definite. In applications 
it often does not depend on coordinates q' . 

Quite independently on the above particular structure, in general equations of 
motion of the system (L, D) have the form: 

Let us quote, although it is not particularly important, that for (|23[) the result- 
ing field equations have the form: 

W _ pr ik F d <l 3 r ik dV r ik H d 1 J . 1 , 
-DV~ eG Fk ^~ G dq~ k ~ k3 ^' a- 1 --,*, (26) 
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where 

Holonomic constraints are described analytically by functionally independent equa- 
tions 

F a (q x ,...,q k ) =0, a=l,...,m, 
and the total system of equations of motion is given by 
D dL dL . , 

ntW-f¥ =Di + Ui > M*V.,**)=0, (28) 
i = 1, . . . , k, a = 1, . . . , m, 

together with the condition that reactions IZi are passive, i.e., do not any work on 
virtual motions compatible with constraints, i.e., such ones that 

(9)^ = (29) 

i 

if equations 

F a (q 1 (t),...,q k (t))=0 ) o = l,...,m, (30) 

are satisfied. Let us notice that q l (t) in (|30)) are quite arbitrary (excepting, of 
course, appropriate smoothness conditions). This implies that reactions are given 
by: 

n(q,v) t = J2^ a (q,v)^i, i = i,...,k. (3i) 

a 

The total system of equations 

DDL M_ D( + EA „^ fc (32) 



Dt dq l dq 1 



F a (q\...,k k ) = 0, a = l, 



..,m, 



consists of (k + m) independent conditions imposed on (k + m) functions of time, 
q l (t), X a (t). And, roughly speaking, A a are to be eliminated, and the general 
solution for the time dependence of q % is to be determined. Although, more pre- 
cisely, in engineering design the quantities A a are very relevant in the analysis of 
the endurance of mechanical factors responsible for constraints (various "rods", 
"threads", etc.). 

This is implicit description. In certain problems the explicit parametric formu- 
lation is more convenient. So, let the constraints W C Q be parametrically given 
by: 

cf = <p i to 1 ,...,y h - m ), i = l,...,k; (33) 
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the parameters /i = l,...,/ = fc — m are independent "proper" coordinates on 
W C Q, i.e., they label "true" degrees of freedom of the constrained system. Then 
our equations of motion may be written as follows, in explicitly irredundant terms: 

D dL w dL w f i r<3A\ 

m^r-^ =Dw ^ M=l,---,/ = fc-m, (34) 

where, obviously, Lyy is the restriction of the original L to the tangent subbundlc 
TW, and D w is the pull-back of D to W C Q, 

L w {y,y):=LL\y),^A, (35) 



2 Nonholonomy and vakonomy 

It is interesting to mention here briefly about non-holonomic constraints, although 
it is not our main subject here [I] -[TO]. 

Non-holonomic constraints based on the natural slide-free rolling on rough sur- 
faces are linear in generalized velocities; in certain situations one deals with "linear 
non- homogeneous" , i.e., affine conditions, 

T a {q, v) = Uai(q)v l + f a (q) = 0, a = l,...,m. (36) 

Of course, holonomic constraints (|30[) may be also formally written in this way, 
then uj a i are given by derivatives of F a -s in (|32[) : 

dF a 

u a i = -Q-j, a = l,...,m. (37) 

There are also mixed situations, when both the motion in Q and instantaneous 
virtual velocities are independently restricted, and one deals with the mixture of 
genuine non-holonomic, and holonomic constraints. Let us do not go here into 
such details, roughly, the true non-holonomy is assumed here to follow from the 
non-integrability of the Pfaff problem: 

w a = Lo ai (q) dq l = , a = l,...,m. (38) 

This means that at any configuration q 6 Q we are given some linear (k — m)- 
dimensional subspace W q £ T q M of virtual velocities. There are various degrees of 
non/integr ability. There are two extreme situations: quasiholonomic constraints, 
when Q is foliated (stratified) by the m-dimensional family of (k — m)-dimensional 
strata (fibers), and the total non-holonomy, when at least locally, at any q G Q 
there is a neighbourhood U C Q such that any of its points may be approached 
from q along a curve compatible with (1361) / (|38[) (compare this with Caratheodory 
formulation of the Second Principle of Thermodynamics) . 
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According to the d'Alembert principle, the reactions maintaining such con- 
straints are ideal, i.e., energetically passive, they do not do any work along virtual 
motions compatible with constraints. The total system of equations of motion is 
given by 

D dL dL „ xn dq i „ 

obviously, the summation convention meant with respect to repeated indices. 

This is again the system of (k + m) equations for (k + m) variables (q l , X a ). But 
in non-holonomic theory some new problems appear which some more than one 
century ago resulted in a big confusion, misunderstanding, simply shame. Namely, 
in the case of variational system, the Lusternik theorem for the variational problem 



5 J L(q(t),q(t))dt = (40) 

constrained by the conditions 

F a (g(t),«(t))=w (ri (g(t))g i (*)=0 (41) 



results in equations: 



D dL dL . dq i „ 



where, however, the reactions are not any longer given by the right-hand sides of 

Lusternik theorem, i.e., confined extremum (confined stationary point) after 
some easy calculations implies that 

a ( dujaj doj a i\ dq j d\x a 



This is something evidently different than the dissipative-free (I3U)) . The obvious 
generalization of (|42p to the dissipative case is 

D dL dL _ „ 

DtW~W = Dl + R - (44) 

Only the second term in (|43l) is analogous to (|3T)|) . when we identify the holonomic 
Lagrange multiplies A a with dfj, a /dt. The system (|43l) & (|44p consists of (fe + to) 
differential equations imposed on (k + m) variables (q l ,[i a ). Let us notice however, 
that this is essentially a system of differential equations, because the Lagrange 
multipliers enter there in a differential, non-algebraic way. In this sense the system 
is more "elastic". And equations of motion (|43l) / (|4"4")l are evidently different than 
(|32|) . Moreover (|43]) / (|44f do not describe rough, slide-free rolling of natural me- 
chanical systems. The " magnetic" -like term controlled by fi is completely strange 
from the point of view of such applications. The system based on flj5D / flUJ) has 
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evidently more degrees of freedom. Such system were called by Russian school 
"vakonomic" (variational axiomatic). In Western literature they are called "sub- 
Riemannian" . Quite unexpected applications were found in financial mathematics. 
Besides, such systems are interesting from the point of view of pure differential ge- 
ometry. It seems that systems of this kind may be used in active control, especially 
when one deals with servomechanisms and computer-aided problems of programme 
motion. Systems based on variational principles and higher-order Lusternik princi- 
ples (higher-order differential, but also integral and functional), have some special 
features interesting from the point of view of energetic balance in active control. In 
particular, this is the case with non-holonomic constraints nonlinear in velocities, 
accelerations and higher-order time derivatives. Certainly they are non-physical 
in natural mechanism of sliding-free constraints, but they are promising from the 
point of view of active control. 

In automatic and active control some first-order constraints nonlinear in veloc- 
ities may be used, e.g., when stabilizing velocity of satellites and space ships. 

If such constraints are given by equations: 

F a (q,v)=0, a=l,...,m, (45) 

then there are good geometric reasons to control the system with the Appell- 
Chetajev reactions: 

a = 1, . . . , m, i = 1, . . . , k. 

In any case, the Appell-Chetajew reactions are geometrically correctly-defined 
(in a manner independent on the choice of coordinates). 

There are also physical reasons to expect the physical utility from the variational 
Lusternik procedure. The corresponding equations of motion have the form: 

§i§-§ = D ' + R - F -^ = °' < 47 » 

with the program forces / reactions given by: 

dF d^dFa d 2 F a d 2 F a dV 

1 M dq i dt 8q l * dq l dql dt * dq l dql dt 2 { ' 

They are also geometrically correctly defined and it is clear that the main term 
has just the Appell-Chetajew form ^j-^|f"- The last term in (|4"5f represents the 
control of inertial properties. 

3 Constraints and symmetries in deformable 
bodies 

What is not clear with constraints if everything is so clear and classical? It seems 
that people doing with continua often do not distinguish between Procedures 1 
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and 2 in Section 1 above. The point is that mechanics of continua is often confused 
with field theory. In spite of using partial differential equations in both disciplines, 
they are something else. There is only one obvious exception. When dealing with 
incompressible (isochoric) media, one usually does not forget about the pressure 
Lagrange multiplier which is just the reaction force responsible for incompressibil- 
ity. But one forgets about the problem in other models, and it is a great merit of 
professor Czeslaw Wozniak [Tl]-[14] that he stressed the problem. Incidentally, the 
author here is very indebted to professor Wozniak for his understanding the prob- 
lem, and his permanent support. First of all, let us notice that from the geometric 
point of view, the most interesting and important constraints are ones implied by 
some symmetry groups. 

These are usually some groups responsible for geometry of the physical space 
or space time, like isometry group, affine group, conformal group, Poincare group, 
Galilei group, etc. Configuration space of various constrained continua very often 
happen to be homogeneous spaces of those groups. One of examples which was 
very interesting for us was afhnely rigid body, i.e., body rigid in the sense of affine 
geometry, homogeneously deformable body. It was the model of internal degrees of 
freedom in Eringen's micromorphic continuum. There are also other interesting ex- 
amples, like molecular vibrations, collective models of nuclei, astrophysical objects, 
geophysical problems, macroscopic elasticity in situations when the wave length is 
comparable with the linear size of the body etc. And here there is plenty of mis- 
understandings. Namely, one often does not distinguish between Procedures 1, 
2 from Section 1. One is often faced with the statement that there is only a small 
family of solutions. The point is that one confuses Procedure 1 with Procedure 
2 and one looks for the special solutions of unconstrained problems, rather then 
on the constrained dynamics with its characteristic reaction forces. An extremely 
strange argument is that both the deformation tensor and the stress tensor are 
constant within the homogeneously deformable body, because of which motion is 
to be trivial. 

Obviously, finite bodies with boundary cannot have constant deformation ten- 
sor, except its interior. And the reaction analysis shows that, in virtue of the 
d'Alcmbert principle, reactions responsible for the affine rigidity do not vanish, 
however, it is their monopole and dipole distributions that vanishes, i.e., therefore 
the total reaction force and the dipole distribution of reactions do vanish. Because 
of this, if the configuration of affine body is given by 

x\r^-t) =r*(t)+^ K (t)a K , (49) 

where r l are coordinates of the centre of mass, Lp l k are internal/relative parameters, 
and a K are material variables, then equations of motion have the form: 

M % = F '- *V^- = "« (50) 
The meaning of symbols is as follows: 
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• M is the total mass of the body, 

M = J dn (51) 

• J KL is the co-moving tensor of inertia in the material space, thus constant, 

J KL = J a K a L d i i{a) (52) 

• center of mass is placed at a K = 0, 

J K = j a K dn(a) = (53) 

• F % is the total force, 



= J T l (a)d^(a 



) (54) 



• N KL is the co- moving dipole offerees distribution, therefore, its spatial/Eule- 
rian components are given by 

N ij = J <p* K tpi L a K a L 'd/i(a) = <p* l j a K a L d^{a). (55) 

Let us quote some alternative forms of equations of motion as balance laws for 
linear momentum p and affine spin K 

dt dt dt dt y ' 

where 

? = M%, KV^K^J™ (57) 
are respectively translational momentum and affine spin. In other words: 

%--«JC», (58) 



n'j = -^f f~ j, (59) 

is an affine velocity, i.e., Eringen's "gyration" 

n A B = ^- 1a Vb^-, (60) 

is its co-moving representation. 
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Let us also quote the following formula: 



^ = W « + 2 f^, (61) 
at ogij 

where the kinetic energy is given by 

T = Ttr+Tmt = T 9 ^nr + 2^—— J ■ (62) 

If Lagrangian is given by 

L = T-V(r i ,<p i K ), (63) 

then 

P^K'j (64) 

are respectively Hamiltonian generators of spatial translations and afhne rotations 
about the centre of mass. 

gij = K H _ K 3i ( 65 ) 
is the spin angular momentum, and 

dS ij 

= N V -N J \ (66) 



dt 



thus, spin is conserved if N lJ is symmetric 

Let us quote a few additional interesting formulas: 

= -p b JbcK ca + F A , ^ = -K AC J CD K DB + N AB , (67) 
dt dt 

where symbol with the capital indices denote co-moving component of physical 
quantities and 

J AC J CB = S A B . (68) 
Another geometrically interesting expressions: 

M— = Mn A B v B + F A , ^RjCA = _n B D n D c J CA + N AB . (69) 
dt dt 

It is clear that we deal here with the system of n 2 + n = n(n + 1) degrees of 
freedom (in the n-dimensional space; physically it is 12, when n — 3) and this is 
just the dimensionality of the general solution, according to the Procedure 2 in 
Section 1. 

There are another interesting problems concerning dynamical afhne invariance 
in mechanics of afhnely rigid bodies and nonholonomic constraints, both usual and 
VAKONOMIC in mechanics of afhne bodies, however, there is no place for them 
here. Some of them are discussed in [2] , [3] , [5] . 
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Let us notice that if the forces are potential, then the dynamical quantities (151)1) 
|57|) are given by 

■dV dV u - 

y dri Y dip k A y ' 

If there exist dissipative forces non-derivable from Lagrangian or Hamiltonian, 
then in addition to (|56|) some additional terms appears. In the simplest case, we 
choose them linear or quadratic in generalized velocities dr l /dt, d^k/dt. 

Let us observe another interesting point, namely, some additional, geometric, 
i.e., group-implied forces imposed onto (j5"6")l . Gyroscopic constraints, or rather 
pseudo-holonomic constraints of rigid motion, consist of the first equation of (|56|) . 



i.e., the assumption that Q a b are respectively g-skew-symmetric and 77-skew- 
symmetric angular velocities in spatial and co-moving representations, 

- -0/ - - 9ik n k i9 ij , h A B = -% A = -vbc^ c dv da , (7i) 

where g is the metric tensor of the physical space and 77 is the material (reference) 
metric. The conditions (|56p are then evidently holonomic and may be written down 
as the conditions of isometry, 

gijip l A(p 3 b = rjAB- (72) 

This is explicitly purely holonomic form. Then the reaction moments Nr are 
evidently symmetric, 

N Rij = N Rji , (73) 

and equations (|56[) are evidently free of explicitly non-specified reactions. Gyro- 
scopic reactions do not vanish, however their full tensor contractions with skew- 
symmetric affine virtual velocities (angular velocities) are vanishing in virtue of 
constraints (|71|) . ()72j) . Taking skew-symmetric part of ()56|) we eliminate reaction 
moments and obtain the effective equations of motion. 

Let us now consider isochoric constraints, i.e., incompressibility. Here one is 
faced with something traditionally very familiar and important in continuum me- 
chanics, first of all in fluids. The traces of affine velocities do vanish then: 

Tr Q = fi% = 0. (74) 

The total contractions of such virtual f2-s with the reaction affine moment Nr 
must vanish: 

NjFSlji = N R i m k l9jk = 0. (75) 
But this means that reactions are pure traces, 

AV, A^\. AY' ,V/-\ A = -TtJVr = - 9ij N R ^. (76) 

n n 
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Therefore, to eliminate the Lagrange multiplier A, we must take the constraints 
condition (|74| (i.e., detip =const) jointly with the g-traceless part of (|56|) itself, 
i.e., explicit, 

One can discuss constraints implied by the linear-conformal group, i.e., gener- 
ated by rotations and dilatations. Then afhne velocity (gyration) has the form: 

fit*,- =u i j + a8 i j, (78) 

where uj % j is the gr-skew-symmetric angular velocity, i.e., it satisfies (|71[) and a is 
an arbitrary real, dilatational parameter, so that 

gij<p i A<P i B = XVAB, A>0. (79) 

Then reaction-free equations of motion consist of the skew-symmetric part of (1561) 
and of the g-trace of that equation, so 

g i ^ i A ^-J AB = g ij N i K (81) 
Reaction moments Nr 13 are symmetric and g-traceless. 

And finally some very interesting example of non-holonomic, but non-VA- 
KONOMIC constraints, when f2 is (/-symmetric. Those are constraints of the purely 
rotation- free motion (the only geometrically correct definition); 

n i j - = - g lk g a tt k i = O. (82) 

Reactions are anti-symmetric, and (|82|) must be joined with the symmetric part of 

(USD, 

^J!£*jAB + ^/^fj AB = + N». (83) 

Think the motion of suspension in a viscous fluid as an example. There are also 
VAKONOMIC models of this kind. But this a quite different story, no place for it 
here. 
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